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(54) Procede et dispositif permettant a un modem de se synchroniser sur un transmetteur de donnees 
numeriques par voie hertzienne en presence de brouilleurs. 

(57) Le procede consiste a estimer (9) la matrice 
de correlation ftxx des signaux recus du trans- 
metteur sur un ensemble de n capteurs de 
reception, a calculer (10) la matrice de correla- 
tion inverse a calculer (11) des vecteurs 
d'intercorrelation f^, entre les signaux recus 
sur ('ensemble des n capteurs et un signal de 
replique d(k) connu, a calculer un critere de 
synchronisation (15) en effectuant les produits 
des vecteurs d'intercorrelation transposes 
conjugues ft«i et de la matrice des vecteurs 
d'intercorrelation fx*, et a comparer (16) la 
valeur du critere obtenu a une valeur de seuil t\ 
determinee pour placer la synchronisation sur 
I'echantillon pour lequel la valeur des criteres 
depasse la valeur de seuil v\. -Pour estimer le 
decalage Doppler, la sequence de synchronisa- 
tion est' decoupee en tranches d'un nombre 
determine de symboles pour effectuer des cal- 
culs de transformation de Fourier rapide sur 
des vecteurs d'intercorrelation determines sur 
chaque tranche. 
Application : Transmission radio de donnees. 
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* 

La presente invention cencerne un proced6 et un dispositif permettant a un modem de se synchroniser 
sur un transmetteur de donnees numeriques par voie hertzienne notamment en presence de brouilleurs. 

Dans les systemes de transmissions de donn6es numeriques par voie hertzienne la synchronisation est 
effectu6e, a partir de sequences de symboles connus du recepteur et intercalees entre les symboles d'infor- 
mation. par conflation entre le signal re$u et un signal de reference construit a partir des symboles connus. Le 
resultat est compare £ un seuil pour decider de la synchronisation. Des systemes permettant la synchronisation 
d'un modem a partir d'un seul capteur et d'un seul recepteur sont par exemple connus de Tarticle D.V. Sarwate et 
M.B. Purseiy intitule "Cross correlation properties of pseudo random and related sequences" paru dans IEEE Vol. 
66, pp. 593-619, May 1980. Mais les techniques mises en oeuvre qui permettent en outre une estimation des ca- 
ract6ristiques du canal de transmission deviennent rapidement inefficaces en presence de brouillage. 

Lorsque le rapport signal a bruit n'est pas suffisant pour la prise de synchronisation par une technique 
mono-capfeur, des techniques visant & accrcltre le rapport signal a bruit et a permettre la synchronisation sont 
utilisees! 

Ces techniques peuvent etre regroup6es eri deux families. Une premiere famille vise a pretraiter le signal 
regu avant d'effectuer la synchronisation du signal pretraite. La seconde famille vise & ameliorer le rapport 
signal & bruit et a effectuer la synchronisation conjointement 

Dans les techniques relevant de la premiere famille il faut egalement distinguer deux categories : 
4 - Une premiere dat^gorie met eh oeuvre un pretraitement mono-capteur de type excision de brouilleurs 
'avant d'effectuer la synchronisation: Ce traitement est connu par exemple de Particle de LB. MILSTEIN 
"Interference rejection Techniques in spread spectrum communications", paru dans Proc. IEEE, Vol: 76, 
pp. 657-661, June 1988. Suivant cet article on cherche a predire la contribution du brouilleur, dont la 
bande est suppos6e plus petite que la bande du signal, sans predire celle du signal utile. Cela permet 
d'ameiiorer les performances de la synchronisation en presence d' interferences bande etroite (type 
FSK). Mais le traitement reste inefficace en presence d'un brouillage occupant toute la bande utilisee 
par le signal utile. Cela est en partie due au fait que les techniques de pretraitement mono-capteur ne 
permettent un traitement des signaux que dans uri espace limits a I'axe des temps et a I'axe des fre- 
quences. 

- La seconde categorie concerne les techniques mettant en oeuvre un pretraitement multi-capteurs de 
rejection, pr6c6dant la prise de synchronisation. L'id6e consiste a mettre en oeuvre un f iltrage spatial, 
discriminant le signal et les brouilleurs par leur direction d'arrivee, le temps, le spectre ou la puissance 
comme cela est d6crit dans Touvrage de RA Monzingo, T.W. Miller, intitule "Introduction to adaptive 
arrays", edite par John Wiley and Sons, New-York, 1980 et la these de P. Chevalier, intitul6e "Antenne 
adaptative : d'une structure lineaire a une structure non Iin6aire de Vol terra". These de doctorat, Uni- 
versite de Paris Sud, Juin 1991 puis de correier le signal de reference avec la sortie de 1'antenne. Les 
techniques de f iltrage spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la puis- 
sance sont generalement sou s-opti males, alors que celles exploitant la connaissance de la direction 
d'arrivee du signal necessitent la mise en oeuvre d'une technique de goniometrie, couteuse et delicate 
dans la gamme des frequences HF notamment du fait des multi-trajets. 
Pour ces raisons, les techniques de la seconde famille; visant a effectuer la rejection de brouilleurs conjoin- 
tement avec la prise de synchronisation, a partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exemple ete propos6e par M. Compton dans son article "An adaptive 
array in a spread spectrum communication system" paru dans Proc. IEEE, vol. 66, pp. 289-298, March 1978, 
pour la synchronisation de signaux PN, a partir d'une architecture anaiogique integrant des "Delay Lock Loop". 

Une seconde technique fut propos6e par L.E. Brennan et l.S. Reed dans leur article intitule "An adaptive 
array signal processing algorithm for communications", paru dans IEEE Trans- Aero-Syst. vol. AES-18, n°1, 
pp. 124-130, June 1982 pour les transmissions numeriques et compietee par D.M. Duglos et R.A. Scholtzdans 
Particle intitule "Acquisition of spread spectrum signals by an adaptive array", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc, vol. ASSP-37, n°8, pp. 1253-1270, August 1989. Ces techniques apparaissent toutefois 
beaucoup plus couteuses que les techniques mono-capteur et sont malheureusement inapplicabies telles quei- 
les en presence des derives Doppler du signal inherentes aux transmissions ionospheriques. 
Le but de I'invention est de pallier les inconv6nients precites. 

A cet effet I'invention a pour objet un proc6d6 permettant & un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistant 

- estimer la matrice de correlation Rxx de signaux re$us sur un ensemble de N capteurs de reception par 
blocs composes d'un nombre d'echantillons d'une sequence de synchronisation transmise par le trans- 
metteur, A 

- a calculer la matrice de correlation inverse R^, 
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- a calculer. des vecteurs d'intercorrelation fxd entre les signaux X(k) recus sur I'ensemble des N capteurs 
. etun signal de replique d(k) cpnnu, j tr . 

-a calculer un critere de synchronisation effectuant ie produit scalaire entre leyecteur d'intercorrelation 
. . rxd et le vecteur obtenu par produit de I'inverse de I a. mat rice de correlation et du vecteur d'inter- 

5 correlation r X d 1 ^ . ; ■ 

- a comparer la valeur du critere obtenu a une valeur de seuil ti determinee pour placer la synchronisation 
sur I'echantillon pour lequel la yaleur des criteres depasse la valeur de seuil T|,.caracterise en ce qu'il 
consiste pourestimer le decalage Doppler inherent auxcommunications entre emetteurs et recepteurs, 
a decouper la sequence de synchronisation en un nombre.determine de M tranches de ko symboles, a 

10 effectuer sur ces tranches des calculs de vecteurs d'intercorrelatioapartielles r xd i (* = 1, M) de telle 

sorte que sur chacune des tranches de k<, symboles, le decalage Doppler fasse tourner la phase d'une 
valeur inferieure a 90°, a effectuer des calculs de transformee de Fourier rapide (FFT) surles vecteurs 
r^i obtenus, et a calculer une valeur du critere de synchronisation a partir des vecteurs r Xd (kAf). obtenus 
en resultat des calculs de transformee de Fourier rapide, et de la matrice de correlation Rxx- 
15 ' ^invention a egalement pour objet un dispositif pour la mise en oeuvre du procede. 

Le procede et le dispositif selon I'invention ont pour avantages qu'ils permettent a la fois d'assurer la syn- 
chronisation tfun modem dans des conditions de^brpyillage severe et en presence de.multi-trajets, d'estimer 
- la-derive en frequence entre emetteur et rfecepteur et ded6tecterles instants d'arriyee de fous les multi-trajets 
associes au signalutile, pourvu que ceux T ci«oientde puissance suff is amment forte par rapport a la puissance 
20 * de bruit de fond. ... . . - . ; 

D'autres caracteristiques'et avantages de I'invention apparaitront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui represented : 

- La figure 1 un schema fqnctionnel d'une antenne adaptatiye selon rinventipn. I 

- La figure 2 un exemple de mise en oeuvre d'un f litre adapte spatial a discrimination par la modulation 
25 pour la realisation d'une antenne adaptative selon Tinvention. 

- La figure 3 un schema d'antenne adaptative pour la mise en oeuvre du procede de synchronisation spa- 
tiale selon I'invention. 

- La f igure 4 un graphe pour illustrer la methode d'estimation de la matrice de correlation de signaux pro- 
venant de N capteurs. 

30 - La f igure 5 un graphe ; pour illustrer le decpupage des symboles de la sequence de synchronisation par 

tranche pour le calcul de vecteur d'intercorrejation.en presence de decalage Doppler. . v 
La figure 6 un graphe montrant revolution de critere de synchronisation mis en oeuvre par I'invention 
en presence d'un decalage Doppler. v 

- La figure 7 un organigramme retracant I'ensemble des etapes necessaires a la mise en oeuvre du pro- 
35 ; cede de synchronisation selon ('invention. t - • . 

^ Le procede de synchronisation mis en oeuvre par Tinvention fait appel, aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer a la realisation d'antennes adaptatives sur replique: 

Sejon la technique de synchronisation mpno-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entre K symboles connus dans une sequence de synchronisation et les K echantillons 
40 qui suiventi'arrivee d'un echantillon n sur le capteur, en normalisant les resultats par la puissance p(n) du signal 
recu. Ces calculs sont effectues en appliquant les relations : 

• y(n) = l2;x(k+n)d*(k) (1) 

: ..• K k=1 . 



t 



• * 



p(n) = ^ i;x(k + n)x*(k + n). (2) ; ' 
K k=i / ; • 

■ c <") =1 Sr (3) ' 

■ ■ . . P(n). 

- x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de replique connu (K symboles places en d6but de trame), 
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- * represente l'op6ration de conjugaisori.sur lesnombres complexes. 

Ces calculs sont r6it6res pour chaque echantillon n, et lorsque le critere C(n) depasse un certain seuil n, 
dependant de ta probability de fausse alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine par le rang d'6chantillon n 0 pour lequel C(rio)> r\. La normalisation par la puissance p(n) permet de 
5 choisir un seuil independant du niveau de signal recu. 

Dans ce proc£d6 la synchronisation est obtenue lorsque le rapport signal d bruit est suff isant Ses per- 
formances dependent en partie du nombre d'6chantilIons pris en compte pour la correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peut etre effect uee. 

L'idee mise en oeuvre dans I'invention consiste alors a effectuer la synchronisation, en sortie d'un traite- 
10 ment d'antenne adaptative, permettant la rejection des brouilleurs. 

Le schema fonctionnel d'une antenne adaptative est represents a la figure 1 , celle-ci comporte un r6seau 
1 de N capteurs 1 1 ... 1 N , uh r6seau de filtres adaptes 2 (constitu6s pardes circuits de pond6rations amplitude- 
phase, des filtres a r6ponse impulsionnelle f inie RIF ou infinie Rll) couple & un processeur de commande adap- 
tatif 3. Le processeur 3 calcule les coefficients \N % a W N du reseau de filtres 2 pour I'adapter au signal utile 
15 et aux brouilleurs, en I'absence de toute connaissance a priori sur les brouilleurs et a partir d'une information 
minimale siir le signal utile. Les coefficients a W N permettent de former un trou dans le diagramme de 
rayonnement du reseau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 
la direction du signal utile. ,.:s* • -r i \>. ' 

Differents types d'information sur le signal utile sont exploits par I'antenne adaptative. Ceux-ci sontd'or- 
20 dre spatial lorsque la direction d'arriv6e du signal utile est connue, d'ordre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en Evasion de frequence, retatifs d la forme 
d'onde lorsqu'une information sur la modulation du signal utile est disponible pour disposer d'une r6plique du 
signal utile en intercalant par exemple des s6quences connues dans la trame utile, d'ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupent une bande plus impor- 
25 tante, bu statistiques pour s6parer les difterentes sources arrivant sur le reseau en exploitant par exemple les 
' statistiques d'ordre 4 des signaux en supposant que les differents signaux sont statistiquement ind6pendants. 
Dans le cas de I'antenne adaptative representee k la figure 2, reformation sur le signalutile est constitute 
par les sequences de symboles connus interca!6es dans la trame entre les sequences de sym boles d'infor- 
mation. Les filtres places en sortie de chacun des capteurs 1! k 1 N sontform6s respectivement par des circuits 
30 de pond6ration amplitude-phase W n & W N des signaux x, (t) ax^t) fournis paries capteurs. Les signaux obtenus 
en sortie des circuits de pond6ration sont appliques 6 un circuit de sommation 4. 

Les caract6ristiques du r6seau de filtres 2 sont calcu!6es en minimisant une Erreur Quadratique Moyenne 
(EQM) estim6e £ entre un signal de r6plique d(k) et le signal y(k) sortant de I'antenne de la fa9on schematise 
& la figure 2 ou les Elements homologues & ceux de la figure 1 sont rep6r6s avec les memes references, 4 est 
35 d6f inie par la relation 

* * * 

K „ K 



40 



i = Z|v(k) - <k)| 2 = Z |w*X(k) - d(k)[ (4) 



k=1 k=1; 



ou: 



- y(k) est la sortie de I'antenne adaptative et est def inie par 

y(k) = W»X(k) (5) 

45 - W est le vecteur poids def inissant I'antenne adaptative, W 6tant le vecteur transpose conjuguS corres- 

pondant. (le signe t en exposant def init Toperateur correspondent.) 

- X(k) est le vecteur form6 par les signaux recus sur les capteurs. 

- d(k) est le signal de r6plique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 a la sortie d'un circuit soustracteur 5. 
so' Le vecteur poids W est calcule par la formule connue de Wiener : 

' W = Rif Xd (6) 



A A 
"XX 

avec : 



. -" - . Rxx = ± ZXCkUHk) (7) 

K K=1 
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matrice de correlation des signaux prove nant des capteurs 1 , a 1 N . (X'(k) etant le vecteur transpose conjugue 
du vecteur X(k)). ' - 

■ ' * - i ■ • ' - ■ • 

. : f xd=^£xwa'(k) 

."■"■'■«.• -* k = 1 <. - ... 

* ■ - 1 - - ^ • ... 

vecteur d'intercorrelation entre X(t) et d(t) (8). . \-- 

L'antenne adaptative ainsi utilisee conduit a une maximisation du rapport signal/(bruit + interference) 
connu sous la designation SNIR abreviation anglo^saxonne de "Signal to Noise +, Interfere nee Ratio" obtenu 
en sortie de Tantenne, avec une possibility de rejection pouvant aller jusqu'a N-1 brpuilleurs, et est appelee 
Filtre Adapte Spatial (FAS). 1 

La recherche de la synchronisation est effectuee selon Tinvention en sortie de I'antenne adaptative en 
discriminant le signal seloh sa modulation, et en utilisant le critere.de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represents a la figure 3 ou les elements homologues a ceux 
des figures i et 2 sont representes avec les memes references. .Le critere de synchronisation idevient en uti- 
lisant les notations precedentes : ' . - * . ■ : 

. ' . . : - C(n) = lf:y(k + n)d(k) (9) ; '. 

■ • . * k=i ... ..; ., ; " 

* * *- * - 1» 1 

' . * * t „ . * : ■ ■ * 

y(k+n) represente la sortie de I'antenne adaptative conduisant a une minimisation de I'eiTeur quadratique.* 
moyenrie ^(n) estimee sur K echantillons entre la replique d(k) et le signal arrivant pris a partir de rechantillqn^ 
n, suivant la relation ^ 



K - - 

l(n) =.. ]T|y(k + n) - d(k)| 2 avec . y(k + n) = W(n)*. X(k + n);. : (id) 

k=1 . i ■ . : . . • ■ - , . 



avec : ' - a A ' - 

W(n) = Rxx- 1 (n).r Xd (n) (11) 
Les differentes correlations sont calculees sur les K echantillons qui suivent I'arrivee d'un echantillon n 

suivant les relations : , 

* ■ 

f ■ ■ . 

• ■ 

* * * * - 

* _ V 

K 



Rxx (n) = ^ Z X(k + n) X* (k + n) (1 2) 



k=1 



K 



?Xd(n)=^i;x(k + n)d*(k) (13) 
* k=l 



Lorsque I'echantillon n correspond a la position de synchronisation, les poids W(n) def inissent le filtre 
adapte, ety(k+n) contientalors le signal utile debarrasse de la contribution des brouilleurs. Dans le cas contrai- 
re, aucun signal n'est correle avec la replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
defend. 

L'expression du critere de synchronisation spatiale se simplif ie en ecrivant : 



K f K \ 

C(n) = -1 XW(n) t X(k + n)d*(k) = WCn) 4 1 XX(k + n) d*{k) (14) 
K k=1 l K k=1 
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C(n) = ^(n)Rxx(n) Jl r X d(n) (15) 
De meme que pour le cas mono-capteur, ce calcul est nbitere pour chaque Schantillon n, et lorsque le critere 
C(n) depasse un certain seuil ry, dependant de la probability de fausse alarme. la synchronisation est decidee 
(la normalisation par la puissance est ici obtenue par rintermediaire de Rxx 1 ). 

L'instant de synchronisation est alors determine par le rang de I'echantiHon ppur lequel C(n 0 ) est superieur 

a une valeur de seuil determinee r\. . J , 

Les calculs correspondant sont effectues sur la figure 3 par le bloc 7 de calcul du critere C(n) et par le 
detecteur de seuil 8, a partir de la sequence de synchronisation d(t) et du signal y(t) correles par le correlateur 
6. 

Cela necessite pour chaque echantillon. d'une part de calculer N correlations entre le signal de rephque 
d(k) et les sig naux x(k) recus par les N capteurs (soit N fois les calculs effectues pour la synchronisation mono- 
capteur). et d'autre part d'estimer puis d'inverserla matrice d'intercorrelation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du filtre adapte 2 sont effectues conjointement. 
15 Lorsque la valeur du critere depasse le seuil calcule. la synchronisation est decidee et le vecteur poids de- 
f inissant le filtre 2 s'en deduit facilement par la formule connue de Wiener. Le vecteur poids W peut ensuite 
etre utilise pour filtrer les echantillons contenant Information. ..... v 

En supppsant que la replique d(t) est de puissance 1 , les correlations tendent vers les valeurs suivantes, 
lorsque le n ombre d'echantiTlons K deviehttres grand : 1 ' : v 1 ' - >- 

20 % a - r~ ~ * ■ ■ • - 

rxd - vies S 



10 



25 



30 



35 



45 



ou: 



- R est la matrice de correlation des interferences et du bruit de fond, 

- S est la matrice du vecteur directeur du signal utile, 

- ic 5 est la puissance du signal utile recu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur C 0 s*6crit encore: ~ : 

i^StR- 1 S = SNIRopt (16 x 
0 = 1 + TisStR-iS 1 + SNIRopt 1 
ou SNIRopt est le SNIR asymptotique obtenu en sortie du filtre spatial adapte 2. 

Lorsque le SNIR est fort, I e critere tend vers t, lorsqu'il est faible, il tend vers 0. ^importance du choix du 
seuil ti, f ixe la probability de detection de la synchronisation. Le seuil depend naturellement de la probability 

de fausse alarme souhaitee. A 

Le calcul du critere de synchronisation C, de f Xd et de Rxx est effectue pour chaque position de la syn- 
chronisation et notamment la matrice Rxx est estimee et invers6e 6galement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est elevee. Mais il est possible de reduire la puissance de calcul 
en estimant la matrice de correlation Rxx moins souvent sur un intervalle de temps plus important (K* 6chan- 
tillons au lieu de K avec K'>K de lafacon representee a la figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisation incluses dans cet intervalle de temps. 
40 Dans ces conditions la matrice Rxx s'ecrit : 



Rxx K . 



k=l k=K n 



+ e (17) 



so ou : 

- B(tfc) est la contribution du bruit et des brouilleurs sur X(t k ), 

- S(tk) est la contribution du signal utile sur Xfo), 

- Kq est l'instant d'apparition du signal utile sur le palier, 

- e est un terme prenant en compte rinterconrelation entre le signal B(t) et le signal utile S(t), qui tend vers 
55 0 lorsque K devient grand. - - - 

Pour la position de synchronisation testee : 
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1 

rXd = K 



K 0 +K. > . . 
£ X(tjd(k) 

k=Kn -■• , 

H • * ' * * I 



(18) 



10 



15 



20 



25 



Dans les conditions asymptotiques (avec K>>1), lorsque le signal utile est synchronise avec la replique 
d(t), Rxx etfxd se simplifient et s'ecrivent : 

K'-Ko + 1 



Rxx = R + 



:r 5 S.S< (19) . . 



A 



K 

^^S..(20) 

Le critere obtenu en estimant la rnatrice de correlation Rxx sur K echantillons pbssede comme valeur 

• .1' * V**.. 

asymptotique, apres-.un calcul utilisant le iemme, diversion matricielle: 



G'n = 



, SNIRopt 



A ^ K' - Ko + 1 SNlRopt 



• (21) 



KV 



Comme 



K '~ Ktt . + 1 < 1. alors C- 0 Src&W ! •*<••*• v : " 



1 1- | *. I • 
*-f t | 4 



1 ■ t " fj. 



La valeur du critere apparaft dans ces cpnditions.augment^e d'un facteur dependant de l* instant d'appa- 
rition du signal utile sur la frequence consideree, ce qui entraine une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul reduite. 

D'autre part, le nombre JC ne doit pas Stre choisi trop grand, sa valeur est determinee pour que sur 
K' echantillons, I'environnement de brouillage puisse etre considere comme stationnaire. 

Le seuil a prendre eh consideration depend de la probability de fausse alarme (Pfa). La (Pfa) est la pro- 
bability qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil n. a lieu en calculant 
Pfa=P(C>T)) lorsque Tantenne recoit uniquement du bruit. 
En presence de bruit, le critere de synchronisation s'ecrit: 
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C = 



V- k=t 



77 2 B k B k 
v K k=1, 



iF £Bk ^ 
/ \ K k=1 



\ 



v 



(22) 
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En suppbsant que B k est un bruit gaussien complexe de puissance a 2 , independant entre les capteurs, 
lorsque K est suffisamment grand, le terme du milieu de la relation (22) devient : 



40 



1 r. 



<\ K 

77 2 B k B k 
V K k=1 



\ 



t 



= a +2 Id 



(23) 



(Id etant la mat rice identity ou unitaire) 

le critere de synchronisation peut alors etre approxime de la facon suivante 



45 



c = 



1 



K 2 a 2 



r K 
£ B k \ 



>y k • ■ 



J 



(24) 
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Soit : 



/.0 



55 



c = 



1 



N 

K 2 rr 2 " ^ 

K a i=1 



k=1 



2>* 



1 



N 



„2 
K CT i=1 



ZKtlf) (25) 
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avec: 



K 

X B kV = Ri + jIi (26) 
k=1 

Bi repr6sente le bruit obtenu sur le capteuri, pour Techantillon k. * " » ■ . 

Rj + jli peut etre approxime par un bruit gaussien complexe de variance K a 2 .-Dans ces conditions R ( et Ij 
sont chacun des bruits gaussiens reels die variance 



En notant T| = R|Ou Ij, C devient: . .v 



2N 



- a. 



C^S(f) (27) 



i/2 2 - 
K a i=1. 

ou T, est une variable aleatoire qui suit une loi normale centree de variance — — . La variable aleatoire Zi 
^ 2 Ti suit done une loi normale centree de variance 1 . 



C devient : 

-j | -2N 

C ^2i<Z Z f (28) 



Comme la suite des Zj peut etre consideree comme une suite de 2N variables aleatoires independantes 
gaussiennes de moyennes nulles et de variances 1 t la variable aleatoire 

r 

2N 

i=1 



suit une loi du % 2 a 2N d eg res de liberte. 

La probability d'avoir C ^ n s'ecrit alors : 



p(C > Tl) = p 



{2N - 



£zf >2Kt! 
Vi=1 J 



= Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils r| obtenues par recherche dans une table donnant 
P(7>A) pour une loi du x 2 a 2N degres de liberte et pour Pfa = 1 Ch 5 . 
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N = 5 




A = 42 




Tl = 0.26 


N = 4 




A = 37 




ti = 0.23 


N = 3 




A.= 33 




ti = 0.21 


N = 2 




A = 28 '< 




71 = 0.18 



Comme entre remission et la reception, le signal transmis peut etre affecte d'un certain decalage en fre- 
quence, ou decalage Doppler, du essentiellemept aux derives en frequences des oscillateurs locaux des 
synthetiseurs & remission et k la reception, ainsi qu'aux pbnditions'de propagation ionosp'herique, il est ne- 
cessaire de compenser ce decalage. Le risque est s'il n'est pas compehse, d'empecher toute prise de syn- 
chronisation. 

En appelant Af 0 , le decalage Doppler entre le signal recu par I'antenne etle signal attendu eten considerant 
la position de synchronisation, le vecteur signal recu peut s'ecrire avec les notations des paragraphes sui- 
vants : 

X(t) = d(t) e^o-t S + B(t) 
La matrice de correlation estimee peut alors s'ecrire : ;j ■ -c -", . * '7 ; 



RXX = ^ Z [V^ d(^)e 2 'i n&S o\ s + B(1k)][^ d^e^^'oVs + B(t k j] (30) 
k=1 



RXX = e + 



^ZkKtkfsS^BO^B^l (31) 
K k=1 L J 



ou £ est une matrice prenant en compte les intercorrelations entre la contribution du signal utile et le bruit b(t). 
Lorsque K»1, dans les conditions asymptotiques, Rxx P©ut s'approximer par : 

Rxx s 7i s '.S S + + a 2 Id (32) 

La relation (32) montre que le decalage Doppler ne modifie pas le calcui de Rxx- En effet, le decalage 
s'elimine naturellement lorsque le signal utile est multiplie par sa quantite conjuguee. 

Le vecteur d'intercorrelation peut s'ecrire : 



F Xd = ^ Zl^T <t*> 2 Afe tk S + B(t k )] d*(t k ) (33) 



p Xd = + £ Zj^Ktk)! 2 e 2j * s] (34) 



OU E 



est ub terme prenant en compte I'intercorrelation entre le bruit B(t) et la replique. Lorsque K»1, r Xd peut 
s'approximer par : 



rxd = 



^Z V^KkTfe 
^ k=1 



' e 



(35) 



La relation (34) montre que le decalage Doppler modifie I'estimation du vecteur d'intercorrelation r Xd et 
que du fait que la phase du terme en exponentiel complexe peut tourner de plusieurs fois 360° sur I'intervalle 
de temps utilise pour le calcui de la correlation, Integration ne peut plus avoir lieu de facon coherente ce qui 
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empeche le critere de synchronisation multicapteurs de converger. Le processus suh/ant permet d'estimer ie 
decalage Doppler et de le compenser pour rendre possible la prise de synchronisation. 

Selon ce processus, les K symboles de la sequence de synchronisation sont decoupes de la fapon repre- 
sentee a ta figure 5 en M tranches de Ko symboles (M Ko = K), de telle sorte que sur chacune des tranches 
de Ko symboles, le decalage Doppler fasse tournerla phase du signal regu parTantenned'une valeur inferieure 
d 90°, af in de permettre une integration coherente de rxd sur chacune des tranches. 

Si r Xdi designe la valeur de rintercorrelation estimee sur la tranche i, celie-ci s'ecrit :. 



r Xdi = 



(i+l)K 0 
^ k=iK 0 +1 



0 2jnAk 

^KkT)! 2 e F * 



S i = 0,...,M-1 (36) 



En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
elementaire consideree, et du decalage Doppler, la relation (36) devient : 



r Xdi = 



^ZV^K(k + iKo)Tf e 

^ k=1 



K, 



2 2jn.^.k 



F e 



e 



• S i = 0 M-1 (37) 



AF, 



Ce qui peut encore s'ecrire : 

' f Xdi = x, e^i" -^ lK o .S i = 0 M-1 (38) 

La relation (38) montre que les vecteurs r x <ji sont des vecteurs, echantillonnes & j^. Leur transformee de 
Fourier discrete est alors : 



M-1 



— 2|7l — — —A 



fxd(Af) = Z W, ? ' Fe/K ° 

1=0 



(39) 



iXd(Af ) =S^e 



M-1 2 j*.^- Af ■ 



Fe - s 



i=0 



(40) 



La relation (40) montre que lorsque Af = Af 0 , le decalage Doppler est compense, ce qui rend maximum la 
valeur du critere de synchronisation r X d(AO t Rxx" 1 rxd(Af). 

Af in d'estimer Af. c'est-&-dire af in d'estimer le decalage Doppler qui conduit a la valeur maximale du cri- 
tere, et af in de reduire le cout de calcul correspondant, une transformee de Fourier rapide est effectuee sur 
les vecteurs rxdi (le nombre de tranches elementaires M etant choisi pourformer une puissance de 2) : ce calcul 
donne : 



M-1 : 

fxdtkAf,)- £ r Xdi e Fe/K o 

i=0 



.ik 



M-1 _ 2 jn.^ 

Z p xdi e M 

i=0 



k = 0 M-1 (41) 



ou Afi est la resolution Doppler telle que : 

M 

Si la resolution est insuf f isante, une amelioration d'un facteur 2 peut etre obtenue en calculant une trans- 
formee de Fourier rapide, sur les vecteurs r x< ii "decalees" en frequence de la fagon suivante : 
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* - . 



kAfj + 



V 



M-1 

Z 

i=0 



rx di e M 



-2^ 

e- M 



I m 
» * * I 



k = 0,... l M-1 (43) 



to 



^operation consiste a multiplier chacun des vecteurs d'intercorrelation r Xd i par le terrhe e-O"™) et a effec- 
tuer ensuite un calcul de transformee rapide sur les vecteurs ainsi obtehus. 

Le critere de synchronisation est ensuite estime sur chacune des cases de calcul des transformees de 
Fourier precedentes par la relation : 
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V 



J 



k = 0,...,2M-1 



(44) 



Bien evidemment, si la resolution est tpujours ( Lnsuffjsante, il^est possible ,de ram6liorer en utilisant le 
meme processus pour calculer des vecteurs. 
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rxd 



kAf-i + 
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30 
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oil n est un entier quelconque. 

Un test sur la case de la transformee correspondant au maximum du critere de synchronisation est effec- 

tue. Si le seuil est depasse pour une valeurko donnee du decalage Doppler, it est alors possible d'aff iner 

mm 

le calcul du decalage Doppler (indispensable pour la demodulation) en calculant la valeur du critere pour dif- 
ferentes valeurs de decalage Doppler entourant ko En effet, a Hssue de cette etape, la resolution sur res- 

timation du decalage Doppler est Pour I'ameliorer, il faut proceder par dichotomie : a chaque iteration et 

de la facon representee a la figure 6, la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux valeurs de decalage Doppler donnant les valeurs maxima du critere lors de l* iteration precedente 
est calcul6e. La recherche par dichotomie est arretee des que la resolution souhaitee sur r estimation du de- 
calage Doppler est obtenue. Apres n iterations de I'algorithme de recherche par dichotomie, la resolution, qui 

etait initiaiement de ^\ est : 

2 



Af i / 2 _ Afi 



(45) 



2 n 2 n+1 

L'enchaTnement des etapes de calcul precedente pour la mise en oeuvre du procede selon invention a 
I'aide de moyens de traitement appropries formes par exemple par des processeurs de traitement du signal 
ou tout dispositif equivalent connu est figure sur I'organigramme de la figure 7. 

II consiste suivant Tetape 9 a estimer la matrice de correlation Rxx suivant la relation (7). Cette matrice 
est estimee sur un bloc coritenant un nombre d'echaritillons K' superieur au nombre d'echantillons de la se- 
quence de, synchronisation. 

L'estimation de Rxx est utilisee ensuite pour toutes les positions de synchronisation incluses dans le bloc. 
Puis la matrice de correlation inverse R^ 1 est cateulee a I'etape 1 0. f xdi qui est def ini par la relation (8) est cal- 
cule a I'etape 11. Pour chaque echantillon correspondant a une position $e synchronisation P dans la suite 
des echantillons .P = 1 a P = K'-K, r x<Ji est calcule pour i = 1 a M sur M blocs de Kq symboles comme defini a 

*■ - * 

la figure 5 suivant la relation 



55 



fxdi = r?- X X(P +K + iK 0 ) d* (K + iK 0 ) 
K '0 K=t 



(46) 
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Des calculs de transformee de Fourier Rapide sont effectues a I' etape suivante 12 sur un nombre deter- 
mine M (M = 8 par exenfiple) de vecteuFs fxdi suivant la relation (41). 

Aux etapes 13 et 14 un calcul de vecteur R^i decale el un calcul de transformee de Fourier rapide suivant 
la relation (43) sont executes pour ameliorer la resolution par un facteur de 2. 

Le critere de synchronisation est calcuie a't'etape 15 en application de la relation (44). 

La comparaison du critere ou seuil t| est effectuee a i'etape 16: Si a I'etape 16 la valeur du critere 




est superieure ou egale au seuil, la position de synchronisation est definie a I'etape 17 sur la position P de 
rechantillon correspondant ^estimation du decalage Doppler peut ensuite etre aff ine a I'etape 18 en calculant 
a nouveau la valeur du critere de synchronisation pour des decalages Doppler entourant le decalage 

precedent. Cette etape peut etre suivie par une etape 19 d'estimation de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation precedemment decrite est effectuee pour chacun des echantillons de 
signal numerise. En pratique, pour ameliorer la precision de la position de la synchronisation il est preferable 
d'effectuer un surechantiltonnage du signal d'un facteur determine L, L etant egal a 4 par exemple. Les trai- 
tements de rorganigramme de la figure 7 et d6crits ci-dessus sont encore utilisables sur le signal ainsi sure- 
chantillonne, mais il est possible d'utiliser une variante dans le procede, afin de reduire la puissance de calcul 
necessaire au traitement. Dans ces conditions, apres avoir calcuie pour chaque echantillon le vecteur rxdi pour 
i = M sur M blocs de Ko symboles suivant la relation (46) uh critere de synchronisation C ( est calcuie sur 
chacun des M blocs par la relation 

... C, = ^RA.f„_ (47) 

: Puis un critere somme de Tensemble des C, est calcuie suivant la relation : 

M 

. . ' ' C^P^Cj (48) 

i=1 

• * ■ • 

Ce critere correspond a une integration non coherente des Ci (cbntrairement a ce qui est fait pour le calcul 
par FFT ou on integre de facon coherente), ce qui se traduit par une remontee du critere sur les echantillons 
de bruit. 

Le calcul du critere a lieu par blocs de L echantillons consecutifs. Sur chacun de ces blocs, rechantillon 
qui donne la valeur maximum de C^P) est selectionne pour effectuer le calcul par transformee de Fourier ra- 
pide necessaire a ('estimation du decalage Doppler comme deja decrit aux etapes 11a 19 de la figure 7. 

Revendicxations 

1. Procede permettant a un modem de se synchroniser sur un transmetteur de signaux numeriques par voie 
hertzienne en presence de brouilleurs consistant a : 

- estimer (9) la matrice.de correlation Rxxde signaux recus sur un ensemble de N capteurs de recep- 
tion (1, ... 1 N ) par blocs composes d'un nombre d'echantillons d'une sequence de synchronisation 
transmise par le transmetteur, A 

- a calculer (10) la matrice de correlation inverse R^, 

- a calculer (11) des vecteurs d'intercorrelation f Xd entre les.signaux X(k) recus sur I'ensemble des N 
capteurs et un signal de replique d(k) connu, 

- a calculer un critere de synchronisation (7) ef fectuant le produit scalaire entre le vecteur d'intercorre- 
lation r Xd et le vecteur obtonu par produit de I'inverse de la matrice de correlation R^ et du vecteur 
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d'intercorrelation rx<j * . . - ■ -, • 

- a comparer (8) la valeur du critere obtenii a une valeur de.seuil Tydeterminee'pour placer la synchro- 
. nisation sur I'echantillon pour-lequel la valeur du critere depasse la valeur de seuil ti, caracterise en 
ce qu'il consiste pour est i mer le decalage bpppler inherent aux communications entre emetteurs et 
recepteurs, a decouper la sequence de synchronisation en uri nombre determine de M tranches de 
ko symbol es, a effectuer (11,1 3) sur ces tranches des caiculs de vecteurs d'intercorrelation partielles 

f Xdi (i = 1 M) de telle sorte que sur chacune des tranches de ko symboles', ie decalage Dopp'er 

fasse tourner la phase d'une valeurinferieure a 90°, a effectuer (12, 14) des caiculs de transformee 
de Fourier rapide (FFT) sur les vecteurs r xd i obtenus, et a calculer une valeur du critere de synchro- 
nisation a partir des vecteurs r Xd (kAf) x obtenus en resultat des caiculs de transformee de Fourier ra- 
pide, et de la matrice de correlation Rxx- 

Procede selon la revendication 1 caracterise en ce que I'estimatibn be la matrice de correlation Rxx est 
obtenu par la relation 

R>Qc=||:xW.x ( T K) 

et le vecteur d'intercorrelation est estime suivant la relation 

- - . . . • ■>»."* ■ . 1 • 

' - r • ' : • ' ■ "if I * • 

. r\ k=1 . 

■ • ■ .... ■ :• ^ . , 

* * - - *•.><_ - 

• " - * f*. * . ". ■ : : ■ * - * . ■ ■ , 

ou p est la position du premier echantillon a partir duquel la correlation' est effectuee entre le vecteur X 
et la replique d. , . . 

Procede selon la revendication 2 caracterise en ce qu'il consiste a choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et a'comparer la valeur maxi- 
mum obtenue a une valeur de seuil determinee pour decider de la synchronisation. 

Procede selon Tune quelconque des revendications 1 a 3 caracterise en ce qu'il consiste pour obtenir une 
meilleure resolution sur le decalage Doppler et une valeur plus importante du maximum du critere; a cal- 
culer (1 3) des vecteurs r xdi decales en frequence et a effectuer un calcul de FFT sur les vecteurs r xdi de- 
cales de Jacon a obtenir la valeur du critere avec une meilleure resolution sur les frequences Doppler. 

Procede selon I'une quelconque des revendications 1 a 4 caracterise eri ce qu'il consiste a surechantil- 
lonner les signaux recus.sur les N capteurs par rapport au rythme symbole, pour ameliorer la resolution 
sur la position de synchronisation.* 

Procede selon I'une quelconque des revendications 1 a 5 caracterisee eh ce qu'il consiste a effectuer le 
calcul de la matrice de correlation Rxx sur un nombre d'echantillons K' superieur au nombre K d'echan- 
tillons de la sequence de synchronisation af in de diminuer la puissance de calcul. 

Procede selon I'une quelconque des revendications 1 a 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir surechantillonn6 d'un facteur L par rapport au rythme symbole, a 
effectuer une etape de pretraitement consistant a calculer pour chaque echahtillon p (p - 1 a L) composant 
un symbole les valeurs C| du critere de synchronisation sur chacune des M tranches de ko symboles et 
a calculer la valeur du critere en effectuant la.somme des criteres Q obtenus pour retenir dans la suite 
du traitement I'echantillon donnant la valeur maximale du critere C = £Cj: 

Dispositif pour la mise eh oeuvre du procede selon Tune quelconque des revendications 1 a 7 du type 
comportant, un reseau de N capteurs (1 ,...1 N ) de signaux de reception couple a un reseau de f ittres adap- 
tes (2) comportant des circuits de ponderations amplitude phase commandes par un processeur de 
commande (3), le processeur de commande (3) etant programme pour 
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estimer la matrice de correlation inverse des signaux re9u X^) et calculer les vecteurs d'inter- 
correlation ?xd entre les signaux regus et le signal de replique d tK) caracterise en ce qu'il comprend 
des mbyens de traitement (5.. .19) pour calculer un critere de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur d'intercorreiation r xd et le vecteur obtenu par le produit de I'inverse 
de matrice de correlation R« et-du vecteur d'intercorreiation ?xd. 

comparer (8) la valeur du critere obtenu a la valeur de seuil t| pour placer la synchronisation sur 
r6chantillon pour lequel ta valeur du critere depasse la valeur du seuil tj, 

et estimer le decalage Doppler inherent aux communications entre emetteurs et recepteurs, en de- 
coupant la sequence de synchronisation en un nombre determine de M tranches de ko symboles, 
en effectuant (11; 13) sur ces tranches des calculs de vecteurs d'intercoifeiation partielles r xd| (i = 
1 , ... t M) de telle sorte que sur chacune des tranches de Kq symboles le d6calage Doppler fasse tour- 
ner k\a phase d'une valeur inf6rieure & 90°, en effectuant (12, 14) des calculs de transform6e de 
Fourier rapide (FFT) sur les vecteurs obtenus, et en calculant une valeur du critere de synchro- 
nisation £ partir des vecteurs r X d(KAfJ, obtenus en r6sultat des calculs de transform6e de Fourier ra- 
pide et de la matrice de correlation Rxx- 
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FIG.1 
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FIG.2 
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Calcul du critere 



C(n)=-^£ y(k + n)d 



(k) 



C (n)= r VJ {n)Ri x (n)rxd 



Xd 
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FIG.3 
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Estimation de R xx sur K' echantillons 
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Estimation de R xx 
sur K 1 echantillons 



Calcul' de R 
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